Elastic Constants of Incommensurate Solid ^He 
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We study the elastic properties of incommensurate solid ^He in the limit of zero temperature. 
Specifically, we calculate the pressure dependence of the five elastic constants (Cn, C12, C13, C33, 
and C44), longitudinal and transversal speeds of sound, and the T = Q Debye temperature of 
incommensurate and commensurate hep ^Ife using the diffusion Monte Carlo method. Our results 
show that under compression the commensurate crystal is globally stiffer than the incommensurate, 
however at pressures close to melting (i.e. P ~ 25 bars) some of the elastic constants accounting for 
strain deformations of the hep basal plane (C12 and C13) are slightly larger in the incommensurate 
solid. Also, we find that upon the introduction of tiny concentrations of point defects the shear 
modulus of ^He (C44) undergoes a small reduction. 

PACS numbers: 67.80.-s,02.70.Ss,67.40.-w 
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I. INTRODUCTION 

An intriguing resemblance between the dependence of 
the shear modulus (SM) and torsional oscillator (TO) 
frequency changes on temperature, amplitude, and con- 
centration of '^He impurities, has been experimentally 
observed in solid ''He at low temperatures ii Crystal de- 
fects are clearly involved in both phenomena, however 
how SM and TO fluctuations are exactly related remains 
yet a puzzle. Day and Beamish identified the stiffen- 
ing of solid helium with decreasing temperature, i.e., in- 
crease of its shear modulus, with the pinning/unpinning 
of dislocations by isotopic impurities. Subsequent exper- 
iments have conflrmed Day and Beamish interpretations^ 
although recent elasticity measurements on ultrapure sin- 
gle crystals seem to suggest that SM variations cannot be 
uniquely understood in terms of mobile dislocations)^ 

Torsional oscillator anomalies were first interpreted as 
the mass decoupling of a certain supersolid fractionj^i^ a 
counterintuitive physical phenomenon that Andreev and 
Lifshitz already proposed in solid helium more than 40 
years ago.'' Supporting this view is the fact that TO 
anomalies appear to occur only in bulk *He.^ Neverthe- 
less, the supersolid interpretation of TO anomalies ap- 
pears to leave open its connection to SM fluctuations 
and diverse theoretical arguments and hypotheses have 
been put forward in an attempt to simultaneously ratio- 
nalize the origins of both anomalies. Anderson, for exam- 
ple, proposes that supersolidity is an intrinsic property 
of bosonic crystals, which is only enhanced by disorder, 
and that the elastic anomaly is due to the generation of 
vortices at temperatures close to the supersolid transi- 
tion.^ From a diametrically opposite standpoint, Reppy 
has argued that the TO behavior is caused by an increase 
of the ■*He shear modulus which mimics mass decoupling 
by stiffening the TO setup. '^'^ Other scenarios somewhat 
more reconciling with the original TO and SM interpreta- 
tions have been also proposed in which for instance mass 
superflow is assumed to occur in the core of dislocations 
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FIG. 1: Shear modulus results obtained for C and IC hep ''He 
expressed as a function of pressure. Experimental data from 
Refs. [29l . Isil . and [s^ are shown for comparison. Solid lines 
represent linear fits to DMC results (see text). 



only when these are static4iiii 

As it can be appreciated, definitive conclusions on the 
roots of SM and TO anomalies remain contentious. In 
a recent paper, Chan et al}3i have shown that for solid 
''He in vycor the nonclassical moment of inertia (NCRI) 
disappears if the TO setup is designed in such a way that 
is completely free from any shear modulus stifl'ening ef- 
fect. This result seems to show that NCRI can be totally 
attributed to elastic effects and not to the existence of 
a supersolid fractionJ^ On the other hand, a recent ex- 
periment in which DC rotation was superposed to both 
TO and SM measures suggested that the cause of both 
anomalies below a critical temperature could have dif- 
ferent microscopic origins. '^^ Also, the source of a small 
peak in the specific heat of ^He^^ at temperatures close 
to that at which TO and SM anomalies appear remains 
yet unexplained. 
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FIG. 2: The shear modulus of IC solid hep He expressed 
as a function of the superfluid fraction at pressures close to 
melting. The dashed line is a guide to the eyes. 



In this work, we study the change in the elastic con- 
stants of solid ■*IIe caused by the presence of small point 
defects concentrations, n^, of 0.5 — 2.0 % . As it has 
been shown, the presence of vacancies induces a finite 
superfluid fraction in the crystal (incommensurate crys- 
tal, IC) so that we can theoretically compare the elastic 
constants of a supersolid with those of the perfect crystal 
(commensurate crystal, C). In particular, we estimate the 
pressure dependence of the elastic constants C^^ 's (Cn, 
C12, Ci3, C33 and C44, where the last one is also known 
as the shear modulus) and derived quantities (the T = 
Debye temperature and transverse/longitudinal speeds of 
sound) of bulk IC and C hep ^He. Our calculations show 
that (i) under moderate and large compressions the C 
phase is globally stiffer than the IC solid, (ii) at pres- 
sures close to melting (i. e. P ~ 25 bars) some of the 
elastic constants accounting for specific strain deforma- 
tions of the hep basal-plane (C12 and C13) are slightly 
larger in the IC crystal, and (iii) the shear modulus dif- 
ference between C and IC ''He crystals is about 10 to 90 
times smaller (in absolute value) than the experimentally 
observed C44 variation caused by the pinning/unpinning 
of dislocations. 

The remainder of this article is organized as follows. In 
the next section, we briefly describe the computational 
methods employed and provide the details of our calcula- 
tions. Next, we present and discuss the results obtained 
and summarize the main conclusions in Sec. IIVI 



II. COMPUTATIONAL METHOD 



FIG. 3: Elastic constant differences between C and IC crystals 
{ps/p = 2 %) of hep *He expressed as a function of pressure. 
The size of the error bars are indicated within the parentheses. 
Positive ACij values indicate softening of the corresponding 
elastic constant in an hypothetical vacancy-induced normal- 
to-supersolid phase transition. 



is usual in DMC, we introduce a guiding wave function 
(gwf) for importance sampling that crucially reduces the 
variance of the statistical estimations. Our gwf model is 
symmetric under the exchange of atoms and correctly re- 
produces the experimental equation of state of solid "^He 
and other quantum crystals ii2r2S We note that DMC en- 
ergies are virtually exact, i.e. are only subjected to sta- 
tistical bias, and ultimately do not depend on the par- 
ticular choice of the guiding wave function. The value of 
all technical parameters, i.e. size of the simulation box, 
population of walkers, and length of the imaginary time- 
step, have been set in order to ensure convergence of the 
total ground-state energy to less than 0.01 K/atom. As 
in previous works, we modeled the *He-^He interactions 
with the Aziz II pairwise potential;22 Further technical 
details of our elastic constant calculations can be found 
in Refs. HHii. 

It is important to stress that DMC Cij estimations es- 
sentially rely on computation of total energies as a func- 
tion of strain thus numerical errors stemming from finite- 
size effects can already be made negligible (i.e. smaller 
than 0.01 K/atom) in computationally affordable simu- 
lation boxes of 24.2 A x 24.2 A x 27.4 A containing 
200 atoms. The IC phase is built by introducing small 
vacancy concentrations of 0.5 — 2.0 % in the crystal. Al- 
though it is well-known that the presence of point defects 
in solid ^He is energetically penalized, ^'^ this route allows 
for simulation of supersolids under tight and controllable 
conditions ,^2,24 ,25 



In this study we employ the diffusion Monte Carlo 
method (DMC), an accurate ground-state approach in 
which the Schrodinger equation of a TV-particle interact- 
ing system is solved stochastically by simulating branch- 
ing and diffusion processes in imaginary timCfi^ As it 



III. RESULTS AND DISCUSSION 

In Fig.[Tl we show the shear modulus of C and IC (with 
Ps/ p = 2 %) hep *IIe expressed as a function of pressure. 
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We find that in both states C44 behaves linearly with 
pressure over all the range of densities considered, i.e. 
0.028 < p< 0.033 A-3 . As one may also see, the shear 
modulus of the C crystal is larger than that of the IC 
solid and the value of the AC44 = C^^ — C44 difference 
increases under compression. It must be noted that the 
numerical uncertainty in our C44 calculations is 5 bar 
thus the predicted AC44 values are rigorously different 
from zero at pressures above 50 bar (see Fig [3]). 

In principle, one may expect that besides pressure 
AC44 variations are also dependent on the imposed frac- 
tion of mass superflow, or conversely, the concentration 
of point defects. However, as we show in Fig [2l such a 
dependence turns out to be rather weak. For instance, 
in the 0<ps/p<3% interval C44 decreases in less than 
the 5 % of its ground-state value and even when an exces- 
sive PsI p value of 7 % is constrained the accompanying 
variation of shear modulus is of just ^ —11 %. Concern- 
ing possible temperature effects, it is well-known that the 
contribution of phonon excitations to the thermal energy 
of solids reduces the speeds of sound by an amount that 
is proportional to T^, so implying a oc dependence 
in the elastic constantsi^ Our zero-temperature conclu- 
sions on AC44, therefore, can be fairly generalized to the 
regime of ultralow temperatures (that is, few mK). In 
fact, the ground-state results reported in this study are 
in very good agreement with those obtained by Pessoa et 
al. for hep ''He at T = 1 K using the path integral Monte 
Carlo Ta.eih.odM^ 

It must be stressed that our C44 results are obtained 
for pure, i.e. with zero concentration of ^He atoms, and 
free-of-dislocations ''He single crystals hence direct com- 
parisons to Day and Beamis h^'^'^^ data obtained in poly- 
crystals turn out to be very complicate. In the light of 
our results, however, one may notice that experimentally 
observed shear modulus variations caused by the pin- 
ning/unpinning of dislocations are of opposite sign and 
about one order of magnitude larger (10 — 20 % in poly- 
crystals and ^ 50 — 90 % in monocrystals^— ) than fluctua- 
tions reported here for hypothetical superfluid mass flows 
of ~ 1 % (see Fig[5]). Consequently, we may conclude that 
if a vacancy-induced normal-to-supersolid phase transi- 
tion occurred in solid helium dislocation-mediated me- 
chanical contributions to C44 would totally overwhelm 
those stemming from mass superflow. Interestingly, Ro- 
jas et al. have recently reported an anomalous softening 
of high quality ultrapure monocrystals in the tempera- 
ture region wherein supersolidity could occur. ^ 

We have also determined the ACij{P) {ij = 
11,12,13,33} deviations describing the response of C 
and IC hep crystals to strain basal plane deformations .~ 
First, we note that all these components also present 
a linear dependence on pressure (see Fig. [31 where nu- 
merical uncertainties are indicated within parentheses). 
Second, ACy slopes are all positive thus implying that 
beyond a certain critical pressure C hep ^He is plainly 
stiffer than the IC crystal. According to our calcula- 
tions, this critical pressure is above 85 bar. Interest- 
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FIG. 4: Top: Calculated longitudinal (L) and transverse (T) 
speeds of sound along the basal plane (b) and c-axis (c) of 
C and IC {ps/p = 2 %) hep ''He as a function of pressure. 
Bottom: Estimated T — Debye Temperature of C and IC 
hep ''He as a function of pressure. Experimental data from 
Refs. and are shown for comparison. 



ingly, C'i2 and C13 are largest, by a small amount, in the 
IC solid at pressures below 50 and 70 bar respectively. 
This outcome comes to show that in an hypothetical low 
pressure normal-to-supersolid phase transition, the final 
supersolid could behave more rigidly than the initial nor- 
mal state under certain strain deformations. Neverthe- 
less we find that the Cee coefficient, which is defined as 
5 ■ {Cii — C12) and can be directly measured in acoustic 
experiments, is always smaller in the IC state. This be- 
havior is analogous to the tendency found for the shear 
modulus although Cee variations are in general larger 
(e.g. at P = 25 bar ACee ~ AC44 whereas at P = 85 bar 
AC66 «2- AC44). 

Finally, Fig. |4] shows the calculated longitudinal and 
transversal speeds of sound {vl and vt) of C and IC hep 
■*He under pressure As one can observe, vl velocities 
along the hep c-axis and basal plane are slightly larger 
in the IC crystal within approximately the same pressure 
range in which AC12 and AC13 deviations are found to be 
negative. Nevertheless, speeds of sound deviations near 
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melting turn out to be so small that in practice these 
could probably not be detected with standard means. 
The same can be concluded about the T = Debye 
temperature for which, as we show in Fig. |31 the cor- 
responding C-to-IC variation is smaller than the typical 
experimental precision. In view of these technical limi- 
tations, it would be very interesting to perform new Cij 
and vl.t measurements on *He at large pressures (i.e. 
P > 60 bar) where C-to-IC differences develop larger. 
To this regard, spectroscopic measurements of the E2g 
phonon mode (i.e. the shear mode corresponding to the 
beating of the two hep sublattices against each other in 
the two orthogonal directions of the basal plane) would 
be particularly helpful since in this type of experiments 
(i) Cij values can be determined with a very small im- 
precision of less than the 2 %, (ii) tiny solid samples are 
needed (i.e. of /im size) thus likely crystal quality issues 
present in SM and TO experiments could be somehow 
alleviated, and (iii) pressure conditions can be efficiently 
tunedi^ 



IV. CONCLUSIONS 

To summarize, we have studied the elastic properties 
of hep solid "^He in a metastable IC state and compared 



them to those obtained for its C ground-state. Our cal- 
culations show that near melting elastic constants Cn 
and Ci2 accounting for specific strain deformations of 
the hep basal plane are slightly larger in the IC crystal. 
At moderate and high pressures, however, the C phase is 
always stiffer than the IC. Also, we find that the appear- 
ance of a finite superfluid fraction (e.g. Ps/ P 1 %) 
caused by the introduction of vacancies unequivocally 
provokes a small decrease of the ^He shear modulus (i.e. 
AC44 ~ 1 %). We argue then that if a vacancy-induced 
normal-to-supersolid phase transition occurred in helium 
crystals containing isotopic impurities and line defects, 
dislocation-mediated contributions to C44 would totally 
overwhelm those stemming from mass superflow. As an 
alternative to usual dynamic experiments focused on the 
search of hypothetical supersolid manifestations, we sug- 
gest to perform spectroscopic measurements of the E2g 
mode of ^He at moderate and high pressures. 
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